Polarization enhancement in two- and three-component ferroelectric superlattices 
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Composition-dependent structural and polar properties of epitaxial short-period CaTiC^/SrTiOs/BaTiOs su- 
perlattices grown on a SrTi03 substrate are investigated with first-principles density-functional theory compu- 
tational techniques. Polarization enhancement with respect to bulk tetragonal BaTi03 is found for two- and 
three-component superlattices with a BaTi03 concentration of more than 30%. Individual BaTi03 layer thick- 
ness is identified as an important factor governing the polarization improvement. In addition, the degree of 
inversion-symmetry breaking in three-component superlattices can be controlled by varying the thicknesses of 
the component layers. The flexibility allowed within this large family of structures makes them highly suitable 
for various applications in modern nano-electro-mechanical devices. 



With modern state-of-the-art epitaxial engineering tech- 
niques, complex perovskite oxide superlattices with three dif- 
ferent consituent layers can be grown to realize a wide range 
of designer materials for a variety of applications^^ Re- 
cently, it was shown that hundreds of atomically thin individ- 
ual layers of CaTi0 3 (CT), SrTi0 3 (ST) and BaTiOg (BT) 
could be grown on a perovskite ST substrate, yielding su- 
perlattices with compositionally abrupt interfaces, atomically 
smooth surfaces and excellent polar properties. 2 ' 3 Since the 
relaxed lattice constants of CT and BT are 0.07 A smaller 
and 0.11 A larger than that of ST (3.905 A), respectively, 
epitaxial strain in such structures results in substantial po- 
larization enhancement. 4 In addition, the inversion symmetry 
breaking present in three-component superlattices 2-3 ' 5 leads to 
even greater flexibility in fine-tuning the ferroelectric proper- 
ties of these materials. 

In this letter, in the spirit of Ref. |4], we investi- 
gate structural and polarization-related properties of two- 
component (or "bicolor") and three-component ("tricolor") 
(ST);(BT) m (CT) n ferroelectric superlattices, where l,n = 
0,1,2 and m = 0, ... ,4, epitaxially matched to a cubic ST 
substrate. We demonstrate that a number of such structures 
have enhanced spontaneous polarization with respect to bulk 
tetragonal BT. We also show that, in agreement with the find- 
ings of Ref. |4], all superlattice layers, including the naturally- 
paraelectric ST ones, are strongly polarized, resulting in a 
smooth polarization profile along the [001] (epitaxial growth) 
direction in the structure. 

All the calculations presented here were performed un- 
der periodic boundary conditions equivalent to the presence 
of short-circuited top and bottom "electrode" layers. A 
plane-wave based DFT-LDA method 6 with ultrasoft pseu- 
dopotentialsi was utilized for the ionic relaxation of the 
1 x 1 x (I + m + n) supercells. 30 Ry wavefunction and 270 
Ry electronic density plane-wave cutoffs were used in all the 
calculations. During the relaxations, the in-plane lattice con- 
stant a of the tetragonal cell was constrained to the theoretical 
lattice constant of cubic ST (3.858 A in this investigation) and 
the out-of-plane lattice constant c was allowed to vary. The 
symmetry in all the calculations was restricted to space group 
P4mm (point group C^), i.e., the ions could move only in 
the [001] direction. The system was considered to be at equi- 
librium when the Hellman-Feynman forces on the ions were 
less than 0.5 x 10~ 3 Ry/bohr and the 0-33 component of the 



stress-tensor was smaller than 0.5 KBar. We used a 6 x 6 x N 
Monkhorst-Pack (MP) mesh 8 for all the Brillouin-zone inte- 
grations, where N — 6/(1 + m + n) for I + m + n < 3, 
N = 2 for I + m + n = 4 and N = 1 for I + m + n > 4. We 
employed the Berry-phase method of the modern polarization 
theory 9 to compute the total polarization of each superlattice. 
The 6 x 6 x 2N MP mesh used in the polarization calcula- 
tions produced well-converged results. 

Calculated structural parameters for the bicolor and tricolor 
(ST);(BT) m (CT)„ superlattices grown on ST, as well as for 
strained bulk CT and BT and unstrained bulk tetragonal BT, 
are presented in Table [I] We find that bulk BT, with a re- 
laxed lattice constant 2% larger than that of cubic ST, expands 
by 5.5% in the [001] direction when pseudomorphically con- 
strained to the ST substrate. Analogously, CT, with relaxed 
lattice constant 1.2% smaller than that of the substrate, con- 
tracts by about 1%. These results are in excellent agreement 
with calculations of Dieguez et al*^ The c/a ratios for the 
bicolor and tricolor superlattices show similar trends corre- 
sponding to BT layers expanding and CT layers slightly con- 
tracting. The tricolors have two sets of tabulated values be- 
cause the lack of inversion symmetry in these systems makes 
ferroelectric displacements along [001] and [001] inequiva- 
lent. It is worth pointing out that the c/a values in the tricolor 
systems are slightly different for the [001] and [001] displace- 
ments, reflecting the presence of a polarization-strain coupling 
in the films. The experimental values for c/a ratios, 11 shown 
in parentheses for some structures, are in excellent agreement 
with the results of our calculations, with the only minor excep- 
tions being for (BT)i(CT)i and tetragonal BT. Experimental 
parameters for the latter were obtained for a thin-film sample 
grown under the same conditions as the superlattice samples. 
The most likely source of discrepancy between the theoretical 
and experimental parameters in both of these cases could be 
the structural defects present in the experimentally grown thin 
films. 

Turning to the polar properties of the superlattices, in 
columns 3 and 4 of Table H] we assemble spontaneous polar- 
izations -P[ooi]> computed with the Berry-phase method, as 
well as polarization enhancement factors with respect to bulk 
tetragonal BT \P\ / Pt-BT for all the superlattices and bulk ma- 
terials studied in this investigation. Starting at the top of the 
Table, we notice that our calculations provide a very high po- 
larization for strained bulk CT. With the aforementioned sym- 
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TABLE I: Structural parameters and polarization in bulk CT and BT, 
and in (ST);(BT) m (CT) n superlattices, pseudomorphically grown on 
ST (with in-plane lattice constant a = 3.858 A). Ppoi] was com- 
puted with the Berry-phase method. \P\/Pt-BT is the polarization 
enhancement factor with respect to the polarization in bulk tetrago- 
nal BT. Experimental values, included for comparison, are shown in 
parentheses where available. The positive direction of the [001] axis 
is from the ST to the neighboring CT layer. 



System 


c/a 


^ooil (C/m 2 ) 


|P|/Pt-BT 


Strained bulk: 








CT 


0.9897 


0.434 




BT 


1.0548 


0.368 




Bicolors: 








(ST)i(CT)i 


1.976 (1.97) 


0.026 (0.000) 


0.11 (0.0) 


(ST)i(BT)i 


2.042 (2.07) 


0.231 (0.059) 


0.95 (0.54) 


(BT)i(CT)i 


2.019 (2.08) 


0.231 (0.085) 


0.95 (0.78) 


(ST) 2 (CT) 2 


3.960 


0.168 


0.69 


(ST) 2 (BT) 2 


4.088 


0.245 


1.01 


(BT) 2 (CT) 2 


4.059 


0.306 


1.26 


Tricolors: 








(ST)!(BT)i(CT)i 


3 ' 018 (3.02) 
3.015 


-°- 200 (0.030) 
0.171 


°- 82 (0.28) 
0.70 


(ST) 2 (BT) 2 (CT) 2 


6 - 049 (6.05) 
6.051 


-0 239 

(0.069) 

0.242 


°- 98 (0.63) 
1.00 


(ST) 2 (BT) 4 (CT) 2 


8.163 
8.165 


-0 295 

(0.133) 

0.298 


L21 (1.22) 
1.23 


(ST) 3 (BT) 3 (CT) 3 


9.083 


-0.260 


1.07 




9.077 


0.244 


1.00 


Bulk tetragonal BT 


1.009 (1.026) 


0.243 (0.109) 





metry restrictions in place, this is the polarization CT would 
have if the zone boundary distortions that are present in its 
actual ground-state crystal structure were suppressed. The 
DFT-based model of Dieguez et a' 10 produces a similar re- 
sult for the strained CT. However, this peculiar property does 
not automatically translate into highly polar CT-containing su- 
perlattices. While CT and BT are both ferroelectric in the 
(constrained) bulk, they behave quite differently in the super- 
lattice geometries. BT is highly polar in monoatomic layer 
systems like (ST)i(BT)i or (ST)i(BT)i(CT)i, and becomes 
even more polar with growing BT-layer thickness m. This 
can be seen from the data for the "constant-concentration" 
(ST) m (BT) m and (BT),„(CT) m bicolor systems, and sub- 
stantiated by the (ST) m (BT) m (CT) m tricolor series as well 
as the BT-rich (ST)2(BT)4(CT)2 system. CT, at variance, 
is practically paraelectric in the monoatomic layer systems, 
with (ST)i(CT)i exhibiting low polarization and (BT)i(CT)i 
showing the same polarization as (ST)i(BT)i. However, as 
the thickness of the CT layer (n) increases, so does its po- 
larization: the (ST)2(CT)2 system is substantially more po- 
lar than (ST)i(CT)i, and (BT) 2 (CT) 2 is more polar than 
(ST) 2 (BT) 2 . In (ST) 3 (BT) 3 (CT) 3 , the ferroelectric displace- 
ments in the CT layer are larger than those in the BT layer. 
These results suggest that CT and BT cells tend to remain 



more polar when assembled into thick layers, as opposed 
to being intermixed with each other and naturally paraelec- 
tric ST. 12 However, in experimentally studied systems, epi- 
taxial growth restrictions place a limit on the values of n 
and to. 13 Above a certain thickness, CT or BT layers re- 
lax to their native in-plane lattice constants and the strain- 
induced polarization enhancement is lost. For example, in the 
(ST)2(BT) m (CT)2 superlattices, the polarization starts to de- 
crease at m ~ 7|2»ii so that a delicate balance between con- 
centration and layer thickness has to be maintained to produce 
superlattices with the strongest polarization enhancement. Fi- 
nally, there is an additional degree of tunability of polar prop- 
erties in the tricolor superlattice family. As shown in the 
bottom part of the Table, the inequivalence of the [001] and 
[001] ferroelectric displacement directions results in two dis- 
tinct values of polarization in these structures. 

Experimental values for the spontaneous polarizations in 
some ferroelectric superlattices, obtained by Lee et a/piil to- 
gether with the associated polarization enhancement factors 
are shown for comparison in Table U The former are consis- 
tently smaller than the polarizations obtained in the calcula- 
tions. However, during the experimental sequences of pol- 
ing and polarization reversal, structural defects and incom- 
plete switching of ferroelectric domains usually do lead to 
substantially reduced values of remanent polarization. Our 
calculations, on the other hand, correspond to an "ideal case" 
(no defects and perfect switching) and thus provide upper- 
bound polarization estimates. Nevertheless, if we compare 
enhancement factors instead of actual polarizations, we see 
certain similar trends in the experimental and theoretical 
data. For example, experimental results for the tricolor sys- 
tems show the same sequence of polarization enhancement 
(from (ST)i(BT)i(CT)i to (ST) 2 (BT) 4 (CT) 2 ) as the calcu- 
lated one. Also the agreement between theoretical and ex- 
perimental polarization-enhancement factors in tricolor sys- 
tems improves with increasing polarization of the system. The 
predicted difference in polarizations induced by inequivalent 
ionic displacements in the tricolor superlattices has not yet 
been observed in the samples of Lee et al due to other sources 
of asymmetry in the experimental geometriespiii 

Since the intricate relation between pseudomorphic strain 
and chemical composition of the layers is the key to the po- 
larization enhancement in the superlattices, we investigated 
the layer-by-layer polarization behavior in these structures. 
Density-functional perturbation theoryii was used to obtain 
Born effective charges for ions in relaxed supercells, and these 
were then used to decompose an "aggregate" supercell polar- 
ization into contributions from individual primitive cells com- 
prising each superlattice. The polarization contribution P\ 
from cell A was estimated using the linearized approximation 

i 9u Xi VX i 

where Z* Xi and Au\i are the effective charge and displacement 
of ion i in cell A, and V\ is the volume of the cell. All polar- 
izations and ionic displacements are constrained to the [001] 
direction, and a superscript zero refers to a non-polar structure 
with ferroelectric displacement removed by "unbuckling" the 
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FIG. 1: Examples of local polarization profiles of the tricolor 
(ST)((BT) m (CT) n superlattices: (a) l,m,n = 1; (b) l,m,n = 2; 
(c) I, n = 2, m = 4; (d) l,m,n — 3. Dashed lines show the abso- 

dBEC- 
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lute value of total average polarization in the supercells (Prooi?)- 



AO and BO planes and moving the "middle" planes back to 
the center of each primitive cell. Equation Q can be used 
with either A or B cation-centered cells. In each case we av- 
erage out contributions of the oxygens in the planes bounding 



elementary cells, while cations at the corners set the volume 
of each cell. The total supercell polarization ^qoi] can ^ e re ~ 
covered from Eq. Q by extending the summation to all the 
ions in the system. Total polarizations obtained by such sum- 
mations are in good agreement with those computed by the 
Berry-phase method. 

In Fig^we present the layer-by-layer polarization profiles 
of the tricolor superlattices from TableU Both A and B cation- 
centered cell profiles are shown for each superlattice, but the 
differences between the two are small. Both types of analysis 
show that all of the layers in the superlattice structures (includ- 
ing the naturally paraelectric ST ones) are strongly polarized, 
and that the local polarization is almost uniform throughout 
the superlattice. As already pointed out in Ref. |4j], the near 
constancy of polarization throughout the superlattice results 
from the minimization of the energy terms associated with the 
polarization charge build-up (V • P) at the interfaces (i.e., it 
is energetically unfavorable for the system to have a large di- 
vergence of the polarization). 

In summary, we have used first-principles methods to study 
structural properties and polarization enhancement in short- 
period perovskite-type thin films containing individual layers 
of CT, ST and BT. We find that the presence of highly po- 
lar CT and BT layers induces substantial polarization in these 
superlattices. The actual degree of polarization enhancement 
strongly depends on an interplay between the concentration 
of ferroelectric components in the superlattice and the pseu- 
domorphic strain it can sustain. Maximum polarization en- 
hancement is achieved by the highest concentration of CT or 
BT assembled into the thickest possible layers that are still 
consistent with a state of full epitaxial strain. Our calculations 
also show that the tricolor superlattices that lack a center of 
inversion have two distinct values of polarization depending 
on the direction of the ferroelectric displacement, which can 
be exploited for fine-tuning of the polar properties. We hope 
that this investigation, by providing additional insight into the 
nature of complex polar perovskite materials, will lead to bet- 
ter understanding and control of their properties in the quest 
for more efficient and environmentally friendly nano-electro- 
mechanical devices. 
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